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Abstract 

We have investigated the correction to the hyperfine structure of heavy mul- 
ticharged ions, which is connected with the nuclear-polarization effect caused 
by the unpaired bound electron. Numerical calculations are performed for 
hydrogenlike ions taking into account the dominant collective nuclear exci- 
tations. The correction defines the ultimate limit of precision in accurate 

theoretical predictions of the hyperfine-structure splittings. 
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During the last decade a significant progress has been achieved in investigations of the 
hypefine structure (hfs) in heavy multicharged ions. To date, accurate measurements of hfs 
sphttings were performed for a number of elements by using optical spectroscopy at the 
experimental storage ring at GSI in Darmstadt and at the Super-EBIT at the Lawrence 
Livermore National Laboratory |I|J^,^,^^,^ . One of the major purposes of these investiga- 
tions consists in testing nontrivial effects of bound-state QED in intense nuclear fields with 
an accuracy on the level of a few percent. Experiments on the hfs splitting allow one also to 
probe internal nuclear structure, in particular, the magnetic moment distribution within the 
nucleus (the so-called Bohr-Weisskopf effect). Though a rough estimate of the latter effect 
can be obtained in the framework of the single-nucleon model, its accurate calculation with 
the use of a microscopic nuclear theory still has to be performed. 

Two attempts to calculate the extended magnetization distribution taking into account 
nuclear many-body corrections should be mentioned here. The first approach is based on 
the dynamic correlation model ||^, which allows to include the excited-core configurations 
in the nuclear wave function. Numerical calculations performed in the framework of this 
model provide accurate values for the nuclear magnetic moments. However, no satisfactory 
agreement with experimental results for the hfs splittings has been achieved ||^. The sec- 
ond approach, which also takes into account the nuclear core polarization by the unpaired 
nucleon, is based on the theory of finite Fermi systems In this case, the correction due 
to the nuclear magnetization distribution calculated for hydrogenlike bismuth turned out to 
be too small compared with the corresponding value deduced from the experiment. 

At present it appears to be unlikely to obtain an accurate theoretical value for the Bohr- 
Weisskopf correction. However, one can eliminate it to a large extend in combined measure- 
ments of the hfs splittings in H- and Li-like ions |]T0[. For the ground state of 209gj80+^ 
independent calculations have lead to similar predictions, 797.1(2) meV [|1T| and 797.15(13) 



meV [jT2[, respectively. While both theoretical results agree with the experimental value 
of 820(26) meV further efforts to measure the predicted splitting with higher precision 
have not been successfull. 



In the present paper we evaluate a correction, which has not been previously considered 
in calculations of the hyperfine structure. The correction is connected with the nuclear- 
polarization effect caused by the unpaired bound electron. While the corresponding contri- 
bution is relatively small in the case of hydrogenlike ions, it turns out to be non-negligible 
for accurate theoretical predictions of the hfs splittings in Li-like heavy ions. Since the 
uncertainty of our calculation is comparable with the magnitude of the nuclear-polarization 
effect itself, the latter cannot be completely eliminated by extracting the Bohr-Weisskopf 
contribution in accordance with Shabaev's idea The effect under consideration sets a 



natural limit up to which one can test bound-state QED, even if a specific difference of the 
hfs splittings of H- and Li-like ions is introduced [|1^ . 



In the following, we shall consider a hydrogenlike ion with a nonzero-spin nucleus, so 
that the total angular momentum of an atom F is defined by the coupling of the nuclear 
spin / with the total angular momentum of an electron j. The hfs energy shift due to the 
magnetic-dipole interaction is given by 

AE^(F) = 1 [F{F + 1) - /(/ + 1) - j(j + 1)] ^n.. (1) 

Here Ann is the hfs constant, which depends on the electron state characterized by the 
standard set of quantum numbers. For electron states with j = 1/2, the hfs splitting 
between the levels with F = I + 1/2 and F = I — 1/2 is just AEnK = (-^ + 1/2)A„k. Since 
the nuclear size is rather small with respect to the radius of the electron orbit, the hyperfine 
structure can be fairly understood in the framework of the external-field approximation. 
This allows to treat the magnetic field of a nucleus as a perturbing potential in calculations 
of the hfs splitting. 

Employing the Dirac equation for the electron in the external field of an infinitely heavy 
pointlike nucleus, one obtains {h = c = 1) p^JT5|JT6|] 



where a = is the fine-structure constant (e > 0), gi is the nuclear g factor, Uj. = n ~ \k,\ is 
the radial quantum number, n is the principal quantum number, n = (j + 1/2) (—1)-'^'+^/^, 



7 = \l f^"^ ~ {(^Zy, N = y {rir + 7)^ + {aZy, and m and nip are the electron and proton 
masses, respectively. Because of various nuclear and QED effects, the experimental value for 
the hfs constant deviates from its Dirac prediction To describe the extended nucleus, 
one usually employs the following parametrization 

AT: = A^,{l-S„,^)il-e^^), (3) 

where the corrections SnK and e„K account for the nuclear charge and magnetic moment 
distributions within the nucleus, respectively |lT7| , |l8| , p!9| , pO| . In addition to the approxima- 



tion (^), the radiative corrections should be taken into account. All leading QED effects 
for the hyperfine structure have been independently calculated by different groups and the 
numerical results are consistent ||21|. For Li- like ions, one can develop a perturbation theory 
with respect to the parameter 1/Z, which accounts for corrections arising from the electron- 
electron interaction [0,^]. The recoil correction, which is due to the finite nuclear mass, is 



negligibly small for heavy ions. 

Here we consider a correction AAnn to the hfs constant (0) due to the nuclear-polarization 
effect, which is caused by the bound electron. More precisely, a core-polarization part of 
the effect is considered only, which is due to collective nuclear excitations. The correspond- 
ing single-nucleon contributions should be generally considered beyond the external-field 
approximation. However, they are assumed to be completely negligible. To describe nu- 
clear polarization, we adopt a relativistic field theoretical approach, which incorporates the 
many-body theory for virtual nuclear excitations within bound-state QED for atomic elec- 
trons p3[ . This approach has been successfully applied in calculations of nuclear-polarization 



effects to the Lamb shift p3| , p^ , p5| and to the bound-electron g factor p6[. To some extent. 



the present formulae are quite similar to those derived in Ref. ||2^. The correction under 
consideration may be represented by the sum of contributions, which are referred to as the 
irreducible, the reducible, and the vertex parts. 

The irreducible contribution to the hfs constant can be written in terms of a multipole 
decomposition as follows: 



2vrj(j + l) 



L>0 



1 2 
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(4) 



where = e/(2mp) is the nuclear magneton, is the one-electron energy, ujl = El — Eq 
are the nuclear excitation energies with respect to the ground-state energy Eq of the nucleus 
and B{EL) = B{EL; — > L) are the corresponding reduced electric transition probabilities. 
The sum over ni runs over the entire Dirac spectrum, while the sum over Ki is restricted to 
those intermediate states, where / + /i + L is even. A two-component radial vector {r\nK,) is 
determined by 



{r\nK,) 



(5) 



where P„K(r) = rgn^{r) and Qn^{r) = rfn^{r), with gnnir) and /^^(r) being the upper and 
lower components of the Dirac wave function, respectively. The radial shape parametrizing 
the nuclear transitions is carried by the functions 



Fr(r) 



An 
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in the case of multipole excitations with L>1 and 
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for monopole excitations, respectively. Here i?o is an average radius assigned to the nucleus 
in its ground state. In Eq. (Q), the matrix element is given by 



a\FL\b) = / drFLir) [Pair)Ph{r) + Qair)Qb{r)] . 



(8) 



The perturbed vector {r\nK), which follows as 



n'^n I I I —21 \ 

{r\nK} = {r\nK), 

„i ^riK ^n'K 



(9) 



can be evaluated analytically by means of the generalized virial relations for the Dirac 
equation |^ (see also Refs. p2|j28[] ). In Eq. (|]), is the Pauli matrix. 
The reducible contribution reads 



AA-^ = E B{EL) E 



L>0 



ni,Ki 



n 2 



(10) 



where A^"^* is the hfs constant given by Eq. (H). The sum I + li + L again should be even. 

The nuclear core-polarization correction to the hfs constant due to the vertex part is 
conveniently represented as the sum of a pole term 
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and of a residual term 



AAl 
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(12) 



respectively. Here AA^"' accounts for the terms with rii = n2 and ki = k,2 in the sums over 
intermediate states. The prime in the sum in Eq. (|r^) indicates that CniKi 7^ ^n2K2 when 
Ki = K2, i.e., the pole contribution is supposed to be omitted. In Eqs. (|TT]) and (|T2]), the 
value I + li + L has to be even. A second condition in Eq. (^) is that the sum li + I2 
should be even as well. The total nuclear core-polarization contribution to the hfs constant 
is determined by the sum of all contributions given by Eqs. (|T0|) , ([Tl|) , and {^). 

In Table |, we present numerical results for some hydrogenlike ions, which are of par- 
ticular experimental interest. The calculations were performed taking into account a finite 
set of dominant collective nuclear excitations. To estimate the nuclear parameters, ul and 
B(EL), in the case of nearly spherical nucleus of ^§361, we employed experimental data 
corresponding to the low-lying vibrational levels in neighbouring even-even isotope of ^gfPb, 
which were deduced from nuclear Coulomb excitation. In the case of giant resonances, we 



utilized phenomenological energy- weighted sum rules which are assumed to be concen- 
trated in single resonant states. In the present calculations, contributions due to monopole, 
dipole, quadrupole, and octupole giant resonances have been taken into account. The infinite 
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summations over the entire Dirac spectrum were performed by the finite basis set method. 



Basis functions are generated via B sphnes including nuclear-size effects ^ . 

Concluding, we have evaluated a correction to the hyperfine structure in heavy mul- 
ticharged ions, which is connected with the nuclear core-polarization effect caused by the 
unpaired bound electron. The correction is exhausted over distances of the order of the 
nuclear size and it is enhanced due to a singular behavior of the hyperfine-interaction op- 
erator. The uncertainty of our calculation can be as large as the nuclear-polarization effect 
itself. It yet determines the natural limitation for testing higher-order QED corrections in 
future experiments aiming for accurate hfs measurements. In the case that the experimen- 
tal value for the ground-state hfs sphtting is used to eliminate the Bohr-Weisskopf effect 
10| , the latter cannot be separated from the effect we have considered. Accordingly, the 
utmost precision for theoretical hfs predictions in lithiumlike heavy ions is determined by 
the nuclear-polarization correction to the ground state in hydrogenlike ions. In particular, 
for the ground state in ^''^Bi^^"'", the nuclear-polarization effect contributes on the level of 
about 0.05 meV. This implies that although numerical calculations of bound-state QED 
corrections provide sufficiently stable results, the conservative estimate of uncertainties for 
the hfs splitting in 209gj80+ quoted in Ref. ||TT| appears to be more realistic rather then the 



one predicted in Ref. |12[. It is also worth noting that because of a similar scaling depen- 
dence of nuclear and radiative corrections upon the principal quantum number, significant 
cancellations of almost all corrections, except for the screened QED contribution, occur in 



a specific difference of the hfs splittings in H- and Li-like heavy ions |T3|. In this case, 
the nuclear-polarization effect might be most relevant for determining the uncertainties of 
accurate theoretical predictions. 
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TABLES 

TABLE I. For various hydrogenlike ions, the nuclear spin/parities 7'^, the nuclear magnetic 
dipole moments n (in units of the nuclear magneton) [29], the nuclear-polarization contributions 
to the hfs constant AAi^ for K-shell electron, and the nuclear-polarization corrections to the 
ground-state hfs splitting AEig are tabulated. Column (a): contributions from low-lying collective 
nuclear modes using experimental values for nuclear excitation energies ujl and electric transition 
strengths B{EL); (b) contributions from giant resonances employing empirical sum rules [30]; (c) 
total effect. The negative value of the nuclear magnetic moment for uranium indicates that the 
level with F = 1 — 1/2 lies above the one with F = 1 + 1/2. The numbers in parentheses are 
powers of ten. 
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H (nm) 






AAis (meV) 






AEu ( 


meV) 


159rp^64+ 3/2+ 


2.014 


0.25(- 


-2)^ 


0.26(- 


-2)b 


0.51(- 




0.10( 


-1) 


165Ho66+ 7/2- 


4.132 


0.38(- 


-3) 


0.27(- 


-2) 


0.31(- 


-2) 


0.12( 


-1) 


175Lu70+ 7/2+ 


2.2327 


0.13(- 


-2) 


0.22(- 


-2) 


0.35(- 


-2) 


0.14( 


-1) 


i«^Re^^+ 5/2+ 


3.2197 


0.25(- 


-2) 


0.66(- 


-2) 


0.91(- 


-2) 


0.27( 


-1) 


203rp|80+ 3^/2+ 


1.62226 


0.16(- 


-2) 


0.29(- 


-1) 


0.31(- 


-1) 


0.31( 


-1) 


209Bj82+ 9/2" 


4.1106 


0.81(- 


-3) 


0.10(- 


-1) 


0.11(- 


-1) 


0.55( 


-1) 


235|y91-^ 7/2" 


-0.38 


-0.36(- 


-2) 


-0.29(- 


-2) 


-0.65(- 


-2) 


0.26( 


-1) 
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